(19) 



3 



Europaisches^Kntamt 
European Patent Office 
Office europeen des brevets 




(12) 



(n) EP 0 718 746 A1 

DEMANDE DE BREVET EUROPEEN 



(43) Date de publication: 

26.06.1996 Bulletin 1996/26 

(21) Numero de depot: 95203453.6 

(22) Date de depot: 12.12.1995 


(51) mtci e G06F 1/035, G06F 7/52, 
G06F 7/544 


(84) Etats contractants d6sign6s: 


(72) Inventeurs: 


DE FR GB 


• Martinez, Georges, c/o Societe Civile S.P.I. D. 




75008 Paris (FR) 


(30) Priorite: 21.12.1994 FR 9415410 


• Chauve, Nicolas, c/o Societe Civile S.P.I.D. 




75008 Paris (FR) 


(71) Demandeurs: 




• LABORATOIRES D'ELECTRONIQUE PHILIPS 


(74) Mandataire: Chaffraix, Jean 


S.A.S. 


Societe Civile S.P.i.D. 


94450 Limeil-Brevannes (FR) 


156, Boulevard Haussmann 


Etats contractants designes: 


75008 Paris (FR) 


FR 




• PHILIPS ELECTRONICS N.V. 




5621 BA Eindhoven (NL) 




Etats contractants designes: 




DE GB 





(54) Multiplieur de booth pour fonctions trigonometriques 



(57) Circuit pour multiplier des donnees conforme- 
ment a un algorithme de Booth (10) dans lequel le co- 
dage (12) des commandes est adapte aux caracteristi- 
ques de symetrie d'une fonction trigonometrique. On 
stocke les valeurs de la fonction relatives a un seul qua- 
drant dans une memoire (1 4) et une unite de conversion 



(11) permet d'utiliser le circuit quel que soit le quadrant. 
Les valeurs de la fonction trigonometrique peuvent §tre 
obtenues par interpolation. Application aux calculs de 
transformees de Fourier ou de transformees DCT. 

Application : transformees de Fourier, transformee 
DCT. 
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Description 

L'invention concerne un circuit pour multiplier des donnees par des valeurs d'une fonction trigonometrique d'une 
variable d'entree, le circuit comprenant : 

des premiers moyens pour executer une multiplication en appliquant un algorithme de Booth, 

des seconds moyens pour gen6rer des commandes appliquees aux premiers moyens, dont des commandes de 

signe, conformement a I'algorithme de Booth, 

des troisiemes moyens pour delivrer les valeurs de la fonction trigonometrique. 



Elle concerne aussi un systeme de traitement numerique utilisant un tel circuit et plus particulierement un systeme 
pour un calcul de transformees de Fourier ou autres. Elle concerne aussi un systeme de transmission numerique 
comprenant un circuit numerique de recuperation de porteuse. 

II est courant d'utiliser des circuits pour le calcul de transformees de Fourier dans le traitement numerique de 
is signal. Une transformed de Fourier directe et une transformee de Fourier reciproque concernent pour la premiere une 
transformation de signaux representes dans le domaine temporel en signaux representes dans le domaine frequenciel 
et pour la seconde une transformation reciproque. Ces transformations qui mettent en oeuvre des fonctions exponen- 
tielles complexes impliquent la manipulation de fonctions trigonometriques. Ces manipulations necessitent de realiser 
un nombre important de multiplications par des fonctions sinus/cosinus. 
20 H est done souhaitable d'utiliser des techniques rapides de multiplication pour diminuer les temps de calcul en y 

associant, si possible, ('utilisation d'un materiel reduit pour limiter les couts de mise en oeuvre. 

Une technique connue repondant a ces specifications est celle de I'algorithme de Booth qui consiste a grouper 
les bits des operandes d'une multiplication par groupes de 3 bits et a generer un nombre reduit de signaux de com- 
mandes, par exemple trois commandes seulement, par groupe de trois bits : 

25 

ZERO : ne rien faire 

SGN commande le signe de ('operation a effectuer 

X2 : determine si le multiplicande doit etre ajoute/retranche (selon SGN) une fois ou 2 fois aux sommes partielles. 

30 Un multiplieur qui met en oeuvre I'algorithme de Booth, a I'aide de ces trois commandes, realise des multiplications 

rapides. 

Neanmoins, des techniques permettant d'ameliorer encore ces performances sont toujours recherchees. 

On connatt le document intitule : "Parallel architecture modified Booth multiplier" de A.R. COOPER, IEE Procee- 
dings Vol. 135, Pt. G., N° 3, June 1988, qui decrit un multiplieur de Booth ameliore. 
35 En regroupant les operations selectionnees par les sorties des encodeurs de Booth et en les faisant operer en 

parallele, il est possible de simplifier le deroulement des operations realisees selon I'algorithme de Booth. Un tel mul- 
tiplieur peut etre utilise dans tout calcul et ne tient pas compte de la particulate de I'application done des donnees qui 
sont a multiplier entre-elles. 

Un objet de l'invention est de simplifier encore la construction d'un multiplieur de Booth en tenant compte de la 
40 specificite des donnees notamment dans I'exploitation des donnees trigonometriques. 

Ce but est atteint avec un circuit pour multiplier des donnees par des valeurs d'une fonction trigonometrique d'une 
variable d'entree, tel qu'il est decrit dans le preambule, dans lequel les troisiemes moyens delivrent des valeurs de la 
fonction trigonometrique sur un domaine limite de la variable d'entree, et en ce que le circuit comporte des moyens 
pour modifier les commandes de signe d'apres des etats logiques de bits de poids fort de la variable d'entree et des 
45 moyens de conversion de la variable d'entree pour faire que le circuit fonctionne avec toute valeur de la variable 
d'entree. 

Ainsi avantageusement, en stockant uniquement les valeurs pour un domaine limite de la variable d'entree (par 
exemple -n/2 <6 <0 pour une fonction cos 0), on reduit tres serieusement la tattle de la memoire. De plus, en faisant 
que la commande de signe soit modifiee apres les moyens qui generent classiquement les commandes, on realise 

50 simplement avec peu de materiel la commande des inversions de signe qu'il est necessaire de faire pour realiser la 
multiplication meme lorsque la valeur 9 d'entree n'est pas dans le quadrant trigonometrique correspondant aux valeurs 
stock6es. En particulier, ceci evite d'avoir a ajouter un multiplieur par - 1 pour effectuer des changements de quadrants 
ce qui consomme a la fois du temps de calcul et de la surface de circuit integre. 

Preferentiellement les troisiemes moyens delivrent des valeurs negatives exprimees en complement a 2. Ainsi 

55 avantageusement, pour un format binaire donne (nombre de bits et repartition des poids binaires), i) est possible de 
representor exactement la valeur - 1 , ce qui ne serait pas possible avec la valeur + 1 , exprimee en complement a deux, 
sans introduire une imprecision de calcul de 1 bit du poids le plus faible. Ceci constitue un avantage important dans 
le cas de multiplications effectuees a I'aide de I'algorithme de Booth decrit 
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Selon un mode particulier de realisation, les moyens qui delivrent les valeurs de la fonction trigonometrique cal- 
culent ces valeurs par interpolation autour d'un nombre limite de valeurs de fonction prealablement memorisees. Ceci 
reduit egalement serieusement la taille de la mernoire. 

Ces differents aspects de I'invention et d'autres encore seront apparents et elucides k partir des modes de reali- 
s sation decrits ci-apres. 

L'invention sera mieux comprise a I'aide des figures suivantes donnees k titre d'exemples non limitatifs qui 
representent : 

Figure 1 : un schema d'un circuit multiplieur selon l'invention. 

Figure 2 : un schema d'un circuit multiplieur selon l'invention muni de moyens pour interpoler des valeurs de fonc- 
io tion entre des valeurs stockees dans un premier bloc-memoire. 

Figure 3 : un graphique d'une partie agrandie d'une fonction F(8) = cos (9). 
Figure 4 : un schema d'un circuit d'interpolation. 

La figure 1 represente un cicuit multiplieur selon l'invention comprenant des moyens 10 MULT pour multiplier un 
multiplicande X par un multiplicateur Y conformement k I'algorithme de Booth, des moyens 12 ENCOD pour generer 
is des commandes conformement k cet algorithme de Booth afin de permettre aux moyens 10 d'operer, et des moyens 
14 MEM pour delivrer les valeurs du multiplicateur Y au multiplieur 10. 

L'algorithme de Booth consiste k ajouter un zero en bit de poids faible puis k regrouper les bits du multiplicateur, 
exprime en complement k 2, en groupes de trois bits consecutifs a partir du bit de poids faible ajoute et k en d6duire 
les operations k effectuer. Deux groupes consecutifs se chevauchent par un bit. Les trois bits d'un groupe imposent 
20 ie type d'operations a effectuer : 
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Code 



(001; 010) 
(000; 111) 
(011) 

(101, 110) 
(100) 



Operation 



- additionner le multiplicande a la somme partielle 

- ne rien faire 

- additionner deux fois le multiplicande k la somme partielle 

- soustraire le multiplicande de la somme partielle 

- soustraire 2 fois le multiplicande de la somme partielle. 
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Ces cinq commandes peuvent etre ramenees a trois commandes seulement : 
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40 



45 



50 



- ZERO - ne rien faire 

- SGN - commande le signe de I'operation a effectuer 

- X 2 - determine si le multiplicande doit etre ajoute/retranche (selon SG N) une fois ou 2 fois de la somme partielle. 

II y a done autant de fois les trois commandes prececentes qu'il y a de triplets de bits dans le multiplicateur Y. Si 
celui-ci contient b bits, le nombre de triplets est egal a I'entier le plus proche soit de b/2 soit de b/2 + 1 selon la parite de b. 

L'invention concerne des modifications de la commande du signe SGN delivree normalement par un encodeur de 
Booth (12). 

Le multiplieur 10 de Booth recoit le multiplicande X et plusieurs lots des trois commandes precedentes. Le multi- 
plicateur Y est code par les moyens de codage 12 ENCOD pour delivrer la pluralite de commandes ZERO, SGN et 
X2. Dans l'invention, le multiplicateur Y represente des valeurs d'une fonction trigonometrique, par exemple cos 6. II 
est obtenu en adressant des moyens de memorisation 14 t par exemple une mernoire MEM, qui ont prealablement ete 
charges par des valeurs de la fonction trigonometrique. L'argument 6 de la fonction trigonometrique est fourni en entree 
du circuit selon un mot binaire B P _ V ...B 0 . En utilisant les proprietes de symetrie de la fonction cos 0, il est possible de 
ne stocker dans la mernoire 1 4 que les valeurs cos 0 correspondant k tJ2 <G <k, e'est-a-dire des valeurs negatives de 
cos 0. Les deux bits de poids fort Bp.-, et B P _ 2 permettent de deduire les valeurs de cos 0 pour les autres valeurs de 0. 

Selon ('invention, on controle la commande de signe SGN de I'encodeur de Booth. Une porte 18 recoit Bp.., et B P _ 2 
et delivre un signal CD. Une famille 1 9 de portes OU - exclusif recoit le signal CD et les signaux SGN issus de chaque 
triplet de bits de Y. En sortie, on obtient de nouveaux signaux de signe SGN' qui sont utilises pour selectionner le type 
d'operation k effectuer dans le multiplieur de Booth 10. Le mot binaire Bp...... B 0 est code en complement k 2. Chaque 

nouvelle commande de signe SGN' est alors deduite de la table I de verite suivante : 
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TABLE I 



Bp-1 


B P _ 2 


CD 


SGN 


SGN' 


0 


1 


0 


0 


0 


0 


0 


1 


0 


1 
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TABLE I (suite) 



15 



20 



1 


0 


0 


1 


1 


1 


1 


1 


1 


0 



Pr6ferentiellement, la memoire 14 stocke des valeurs negatives de cos 6 done avec k/2 <6 <7t ou 7t<8 <3n/2. 

Lorsque la valeur 0 qui arrive en entree du circuit n'est pas dans le quadrant correspondant aux valeurs stockees 
dans la memoire, alors une unite de conversion 11 transforme la valeur 0 d'entree d'apres les etats logiques des bits 
de poids 6lev6 Bp. v B P . 2 . Selon le contenu de la memoire 14 (sinus ou cosinus), la conversion va consister a convertir 
io 0 en 7i/2 - 0 et a modifier eventuellement le signe selon le quadrant concerne. 

Pour reduire encore les moyens materiels, les moyens 14 qui delivrent la fonction trigonometrique peuvent operer 
par interpolation autour d'un nombre limite de valeurs de la fonction trigonometrique prealablement memorisees (figure 
2). 

Le mot binaire B P _ V .. B 0 est alors separe en trois champs : 



un premier champ I forme par des bits de poids faible du mot binaire, 
un second champ II forme par des bits medians du mot binaire, 

un troisieme champ III, forme par les bits B P _ V .. B P _ 2 , qui sert a controler les commandes de signe comme expose 
auparavant. 



Les premier et second champs I et II servent a adresser (figure 2) une unite de conversion 11 CONV qui transforme 
les valeurs d'entree 0 en valeur de 0 situee dans le quadrant pour lequel on a memorise prealablement un nombre 
Iimit6 de valeurs de la fonction trigonometrique dans un premier bloc-memoire 142. Les valeurs stockees peuvent etre 
prealablement corrigees par des coefficients correctifs specifiques destines a reduire des erreurs engendrees par la 

25 methode d'interpolation utilisee. 

L'unite de conversion adresse aussi un second bloc-memoire 144 dans lequel on a memorise prealablement des 
coefficients d'interpolation calcules pour les valeurs de la fonction trigonometrique qui ont ete stockees dans le premier 
bloc-memoire. Les coefficients d'interpolation peuvent etre une derivee vraie ou une derivee approximative formee par 
exemple par la pente moyenne du segment joignant deux valeurs stockees consecutives de la fonction. 

30 Ainsi en adressant les deux blocs-memoires avec le champ median II (avec ou sans conversion selon le cas), par 

exemple on peut lire la valeur de la fonction cos 0 (avec corrections) et celle des coefficients d'interpolation pour cette 
valeur 0. 

Le champ I de bits de poids faibles est alors exploite pour calculer une quantite interpolee liee a I'increment fourni 
par le champ I. 

35 La figure 3 represente une partie d'une fonction F(8) = cos(0) connue par un nombre limite de valeurs F(0 nk ) = cos 

(0 nk ) [point A], F(0 (n+1)k ) - cos(0 ( ,^ 1)/f ) [point B], pour des valeurs 0 nk , 0 (n+1 ) k de la variable 0. Sur le graphique n est 
I'indice courant des points A, B..., avec 0 <n<-£-1 , k etant un nombre entier determinant le facteur d'interpolation, avec 

Lorsqu'une valeur quelconque 0 nk+j de la variable est presentee en entree du circuit, celui-ci calcule, par interpo- 
40 lation, une valeur interpolee proche de la valeur F(0n fr+/ ) [point C]. Sous cette forme, I'indice i est un indice courant 
variant de 0 a k-1 afin de diff^rencier k valeurs intermediates entre 0 nk et 0 (n+1 ) k . 

Pour repr6senter N valeurs de la variable 0, il faut log 2 N bits. Une valeur quelconque 0 nk+j de la variable est 
representee par un mot binaire pour lequel on ignore les deux bits de poids forts du debut pour ne considerer que les 
valeurs appartenant a un seul quadrant. Le mot binaire 1 3 (figure 2) comprend un premier champ 1 3A comportant des 
4S bits de poids fort et un second champ 1 3B comportant des bits de poids faible. 

Selon I'invention, Interpolation est effectuee en exploitant un nombre N/k de valeurs connues de la fonction. Ces 
N/k valeurs connues sont pr6f6rentiellement stockees prealablement dans un bloc-memoire 142. Les valeurs connues 
F(0 nk ) sont celles qui correspondent aux log 2 N/k bits de poids fort de la variable d'entree. 

Pour les N/k valeurs connues de la fonction (points A, B, etc), selon I'invention, on stocke egalement, preferen- 
so tiellement dans des moyens de memorisation 144, N/k coefficients d'interpolation C(0 nk )F Ceux-ci peuvent etre des 
valeurs de derivee F'(0 nk ) de la fonction F(0 nk ). Ces valeurs de derivee peuvent etre des valeurs moyennes de derivee 
calculees autour de chaque point connu A, B, etc. Preferentiellement, ces valeurs de coefficient d'interpolation sont 
egales a la pente des segments reliant deux points connus consecutifs. 

Ces valeurs de coefficient d'interpolation peuvent egalement etre formees par des valeurs vraies de derived cal- 
55 culees aux points connus A, B, etc. 

Pour limiter les imprecisions, on stocke dans le bloc-memoire 142 des valeurs connues de la fonction prealable- 
ment corrigees par des termes correctifs. Le bloc-memoire stocke done des valeurs F(0 nk ) + cor(0 nk ). Le terme cor 
(0 nk ) est prealablement calcule pour chaque point connu pour minimiser Cerreur quadratique moyenne entre les valeurs 
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interpolees et la fonction trigonom6trique parfaite. 

Lorsqu'une valeur d'entr6e 8 nk+j est presentee a I'entr^e du circuit, le champ 13A de bits de poids fort adresse 
(connexion 15) le bloc-memoire 142 et les moyens de memorisation 144 : 

s - le bloc-memoire d6livre F(e nk ) + cor(0 nk ) constituant la valeur approchee de la valeur interpose a determiner 
(connexion 17), 

les moyens de memorisation delivrent les coefficients d'interpolation C'(6 nk ) (connexion 16). 

Considerons le cas ou la fonction F(6 nk ) est une fonction cosinus et ou les valeurs de coefficients d'interpolation 
10 sont deduites de Ea pente d'un segment y = a x+b. Ce segment passe au voisinage des deux points connus A et B 
pour indiquer que des coefficients de corrections interviennent dans le calcul des valeurs interposes. 
Une valeur approchee y de cos (6 nk+j ) [point D] est alors telle que : 

cos e n . k -cos e nk 

y= 9 , e « W <W + *nk 

15 ^(nk+k) ^nk 

avec 

\,= 2 -R (nk+i) 

20 0< /< /c-1 0< n < ^ -1 

Apres des manipulations trigonometriques connues, on obtient : 

y = f- 2 sin <2"+l>*tt , sin Ml * + cos nk 

[ N N J k N 

25 Pour obtenir une valeur interpolee cos (9 nk+ j) plus proche de la valeur exacte cos (6 nk+ j), on introduit un terme de 

correction cor (8 nk ). Cette valeur interpolee est situee entre les ordonnees des points D et C. Chaque terme de cor- 
rection cor (9 nk ) est determine en rendant minimum la somme X(c6s B nk+i - cos 8 n ^ +/ ) 2 pour 0 < i < k-1 . 
La valeur approchee cos (G nfc+ /) s'ecrit alors : 
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c6s (6^,) = {-| sin (1^1 .* n ). sin J/jKj J d * cos (_2rt nk j , cor ^ 



Le circuit comprend des moyens de calcul 145, par exemple un multiplieur/accumulateur qui recoit C'(6 nk ) et le 
35 champ 13B forme des bits de poids faible de la valeur d'entree 6 nk+i . 

Dans une premiere etape, les moyens de calcul 145 calculent le resultat partiel : 

c(e nk )i 

et dans une s^econde etape, ils ajoutent ce resultat partiel a la valeur approchee F(9 nk ) + cor(8 nk ) pour delivrer la valeur 
40 interpolee : F(e n/f+/ ) soit : 

F(e ofc+/ ) = [(C'te^) . ,1 + F{Q nk ) + cor(Q nk ) 

Ainsi, il est possible de reduire la taille de la seconde memoire en ne stockant qu'un nombre limite de valeurs sur 
un nombre reduit de bits tout en obtenant une precision elevee sur le resultat, ce qui est particulierement interessant 
45 dans la realisation de circuits integres. 

Dans le cas ou ils delivrent des valeurs interpolees, les troisiemes moyens comprennent : 

le premier bloc-memoire 142 pour m§moriser un nombre limite de valeurs negatives de la fonction correspondant 
au second champ de bits, a des adresses issues du second champ de bits, 
so - le second bloc-memoire 144 pour memoriser des coefficients d'interpolation associ6s aux valeurs memorisees 
dans le premier bloc-memoire, a des adresses issues du second champ de bits, 
des moyens d'interpolation 

pour lire, dans le premier bloc-memoire, une valeur approchee de la fonction en I'adressant par le second 
55 champ de bits de la variable d'entree, 

pour lire, dans le second bloc-memoire, le coefficient d'interpolation en I'adressant par le second champ de 
bits de la variable d'entr6e, 

pour calculer une quantite interpolee en multipliant le premier champ de bits de la variable d'entr£e par le 
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coefficient d'interpolation adresse, la quantite interpolee etant totalisee avec la valeur approchee pour donner 
la valeur interpolee. 

Le circuit peut etre utilise* dans un systeme de traitement numerique de donnees pour effectuer des calculs de 
5 transformees de Fourier, de transformees en cosinus discretes, ou d'autres transformees faisant appel a des fonctions 
trigonometriques. 

Le circuit peut ausst etre utilise dans un systeme de transmission numerique dans lequel le circuit selon I'invention 
entre dans la realisation d'un circuit numerique de recuperation de porteuse. Dans ce cas, la variable d'entree est une 
erreur de phase G k partir de laquelle le circuit multiplie des valeurs cos 0 et sin 6 par un multiplicande forme par un 
10 symbole qui doit subir une correction de phase pour synchroniser la recuperation de porteuse. 

Revendications 

15 1 . Circuit pour multiplier des donnees par des valeurs d'une fonction trigonometrique d'une variable d'entree, le circuit 
comprenant : 

des premiers moyens pour exScuter une multiplication en appliquant un algorithme de Booth, 
des seconds moyens pour generer des commandes appliquees aux premiers moyens, dont des commandes 
20 de signe, conformement k I'algorithme de Booth, 

des troisiemes moyens pour delivrer les valeurs de la fonction trigonometrique, 

caracterise en ce que les troisiemes moyens delivrent des valeurs de la fonction trigonometrique sur un domaine 
limite de la variable d'entree, et en ce que le circuit comporte des moyens pour modifier les commandes de signe 
25 d'apres des etats logiques de bits de poids fort de la variable d'entree et des moyens de conversion de la variable 

d'entree pour faire que le circuit fonctionne avec toute valeur de la variable d'entree. 

2. Circuit selon la revendication 1 caracterise en ce que les troisiemes moyens deiivrent des valeurs negatives de 
la fonction trigonometrique. 

30 

3. Circuit selon les revendications 1 ou 2 caracterise en ce que les commandes de signe sont inversees lorsque les 
deux bits de poids fort de la variable d'entree sont identiques. 

4. Circuit selon la revendication 3 caracterise en ce que les troisiemes moyens calculent la valeur de la fonction 
35 trigonometrique par interpolation autour d'un nombre limite de valeurs de fonction prealablement memorisees. 

5. Circuit selon la revendication 4 caracterise en ce que la variable d'entree etant representee par un premier champ 
de bits de poids faible, par un second champ de bits de poids intermediate, et par un troisieme champ de bits de 
poids fort, les troisiemes moyens calculent, par interpolation, une valeur interpolee de la fonction trigonometrique 

40 pour la valeur de la variable d'entree, les troisiemes moyens comprenant : 

un premier bloc-memoire pour memoriser un nombre limite de valeurs negatives de la fonction correspondant 
au second champ de bits, k des adresses issues du second champ de bits, 

un second bloc-memoire pour memoriser des coefficients d'interpolation associes aux valeurs memorisees 
45 dans le premier bloc-memoire, k des adresses issues du second champ de bits, 

des moyens d'interpolation 

pour tire, dans le premier bloc-memoire, une valeur approchee de la fonction en I'adressant par le second 
champ de bits de la variable d'entree, 
50 • pour lire, dans le second bloc-memoire, le coefficient d'interpolation en I'adressant par le second champ 

de bits de la variable d'entree, 

pour calculer une quantite interpolee en multipliant le premier champ de bits de la variable d'entree par 
le coefficient d'interpolation adresse, la quantite interpolee etant totalisee avec la valeur approchee pour 
donner la valeur interpolee. 



Systeme de traitement numerique de donnees caracterise en ce qu'it comprend au moins un circuit selon une des 
revendications 1 & 5 pour un calcul de transformee de Fourier ou de transformee en cosinus discrete ou de trans- 
formee faisant appel k une fonction trigonometrique. 
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Systeme de transmission numerique caracterise en ce qu'it comprend au moins un circuit selon une des revendi- 
cations 1 a 5 pour realiser un circuit numerique de recuperation de porteuse. 
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